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Abstract 

A  capacity  fade  prediction  model  has  been  developed  for  Li-ion  cells  based  on  a  semi-empirical  approach.  Correlations  for  variation 
of  capacity  fade  parameters  with  cycling  were  obtained  with  two  different  approaches.  The  first  approach  takes  into  account  only  the 
active  material  loss,  while  the  second  approach  includes  rate  capability  losses  too.  Both  methods  use  correlations  for  variation  of  the  film 
resistance  with  cycling.  The  state  of  charge  (SOC)  of  the  limiting  electrode  accounts  for  the  active  material  loss.  The  diffusion  coefficient 
of  the  limiting  electrode  was  the  parameter  to  account  for  the  rate  capability  losses  during  cycling. 
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1.  Introduction 

Li-ion  batteries  have  replaced  Ni-MH  and  Ni-Cd  systems 
for  several  applications  that  require  high  power  densities. 
Numerous  studies  have  reported  the  capacity  fade  of  these 
batteries  under  extended  cycling.  For  aerospace  and  other 
advanced  applications,  it  is  essential  to  quantify  the  capac¬ 
ity  loss  for  a  given  cycling  protocol.  Real-time  cycle  life 
testing  becomes  prohibitively  expensive  if  one  considers  the 
number  of  different  variables  that  can  change  such  as  depth 
of  discharge  (DOD),  charge  and  discharge  rates  and  cycling 
temperature.  Thus,  accelerated  cycle  life  testing  and  devel¬ 
oping  correlations  based  on  this  data  becomes  very  critical 
for  estimation  of  the  capacity  fade  of  the  cell. 

Bloom  et  al.  [1]  reported  accelerated  calendar  and  cycle 
life  study  of  18,650  cells.  Useful  life  according  to  these 
authors  was  strongly  affected  by  temperature,  time,  state  of 
charge  (SOC)  and  change  in  state-of-charge  (ASOC).  To 
estimate  the  percentage  of  power  loss  the  data  were  fit  with 
the  general  equation: 

Q  =  A  exp  (!) 

where  Q  represents  the  percentage  power  loss,  t  (weeks)  the 
time,  T  (K)  the  temperature,  A  the  pre-exponential  factor, 
£a  (J)  the  activation  energy,  and  z  the  adjustable  parameter. 
Similarly,  Inoue  et  al.  [2]  used  the  calendar  life  data  to  pre¬ 


*  Corresponding  author.  Tel.:  +1-803-777-7314;  fax:  +1-803-777-8265. 
E-mail  address:  popov@engr.sc.edu  (B.N.  Popov). 


diet  the  cycle  life  of  the  battery.  According  to  these  authors 
the  corrected  cycling  capacity  loss  data  correlated  linearly 
with  the  cell  internal  resistance.  The  above  models  are 
empirical  in  nature  and  do  not  present  any  insight  into  the 
processes  that  contribute  to  the  capacity  and  rate  capability 
losses  during  cycling. 

The  capacity  fade  and  prediction  of  the  battery  discharge 
characteristics  with  cycling  could  be  obtained  by  using  a 
semi-empirical  approach  or  by  developing  a  first  principles 
model.  A  survey  of  the  literature  indicates  that  for  other  bat¬ 
tery  systems  capacity  fade  analysis  has  relied  on  developing 
empirical  correlations  that  are  specific  to  the  experimental 
conditions.  Hafen  and  Corbett  [3]  developed  an  empirical 
correlation  for  cycle  life  prediction  of  Ni-Cd  cells  based  on 
the  experimental  cycling  data.  The  cycle  of  Mth  failure  has 
been  derived  as 


Cycle  |  Mth  failure 


/  M\  0,229  1  545 

9.71  (  —  j  (DOD)-1'545 

1  3843 


xexp 


T(K) 


where  N  is  the  number  of  cells  in  battery  and  DOD  is  the 
depth  of  discharge. 

Bro  and  Levy  [4]  developed  an  Arrhenius  type  relation 
to  explain  the  capacity  loss  of  Li-SOCl2  cells  with  storage 
time  and  aging  temperature,  which  is  expressed  as 

C(t,  T)  =  Co  -  k0te~(mT)  (3) 


where  C(t,  T)  is  the  capacity  loss  at  time  Y  and  tempera¬ 
ture  T  and  the  capacity  loss  at  t  =  0  is  defined  as  Co  =  0. 
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Nomenclature 

a 

specific  surface  area  of  porous 
electrode  (m2/m3) 

c 

concentration  of  Li  or  Li+  ions  (mol/m3) 

D 

diffusion  coefficient  (m2/s) 

E 

Faraday’s  constant,  96,487  (C/mol) 

io 

exchange-current  density  (A/m2) 

Lpp 

applied  current  density  (16.54  A/m2) 

J 

local  volumetric  current  density  (A/m3) 

k 

rate  constant  of  electrochemical 
reaction  ((A/m2)/(mol/m3)3y/2) 

L 

length  of  the  cell  (m) 

MW 

molecular  weight  (mol/kg) 

N 

cycle  number 

Q 

capacity  (mAh) 

r 

radial  coordinate  (m) 

R 

particle  radius  (pan) 

Rf 

film  resistance  at  the  electrode-electrolyte 
interface  (Q  m2) 

universal  gas  constant,  8.314  (J/mol) 

t 

time  (s) 

t+ 

transference  number  of  Li+  ions  in  the 
electrolyte 

T 

temperature  (K) 

U 

local  equilibrium  potential  (V) 

X 

coordinate  across  the  cell  thickness  (m) 

Greek  letters 

anodic  and  cathodic  transfer  coefficients 
of  electrochemical  reaction 

£ 

volume  fraction  of  a  phase 

0 

local  potential  of  a  phase  (V) 

1 1 

local  over-potential  driving 
electrochemical  reaction  (V) 

K 

conductivity  of  electrolyte  (S/m) 

0 

state-of-charge 

a 

conductivity  of  electrode  (S/m) 

P 

density  of  active  material  (kg/m3) 

Subscripts 

1 

solid  phase 

2 

solution  phase 

n 

negative 

N 

cycle  number 

P 

positive 

s 

separator 

th 

theoretical 

Superscripts 

0 

initial 

eff 

effective 

max 

theoretical  maximum 

ref 

reference 

The  parameters  namely  acceleration  factor,  ‘ko’a nd  activa¬ 
tion  energy  ‘E  were  obtained  from  the  least  squares  fit  from 
a  plot  of  the  logarithm  of  the  capacity  loss  obtained  experi¬ 
mentally  versus  inverse  of  the  aging  temperature. 

Researchers  at  SAFT  [5]  developed  a  life  duration  model 
for  the  Ni-MH  cells  that  involves  an  empirical  correlation 
for  rate  of  corrosion  given  as 


d 


corrosion 


(t,  T,  %) 


dt 


de  ds 


dt  ’  dt 


,  e,  S,  E 


a 


where  t  is  the  time,  T  the  cell  temperature,  percentage  de¬ 
notes  the  discharged  capacity,  S  the  initial  alloy  surface 
area,  Ea  the  activation  energy  and  e  the  thickness  of  the 
corrosion  layer.  In  this  paper,  we  plan  to  adopt  a  similar  ap¬ 
proach  to  derive  semi-empirical  correlations  for  the  capacity 
fade  of  Li-ion  battery  systems.  The  main  drawback  of  this 
methodology  is  its  system  specific  nature.  In  other  words, 
the  semi-empirical  correlations  developed  for  one  type  of 
Li-ion  cell  of  specific  geometry  and  electrode  material  need 
not  be  the  same  for  other  type.  Even  though  the  parameters 
considered  for  capacity  loss  may  be  common  for  all  Li-ion 
systems,  the  trend  in  the  variation  of  these  parameters  with 
continuous  cycling  may  not  be  the  same.  Thus,  in  order  to 
develop  a  mathematical  model  for  capacity  fade,  it  is  critical 
to  understand  the  various  mechanisms  involved  in  capacity 
loss.  Our  previous  research  on  Li-ion  systems  of  different 
geometries,  cathode  materials  and  varying  cycling  condi¬ 
tions  resulted  in  a  better  understanding  of  the  capacity  fade 
mechanisms  involved  in  each  of  the  systems  studied  [6,7]. 

This  paper  primarily  focuses  on  developing  semi-empirical 
correlations  to  account  for  the  capacity  fade  in  Li-ion  cells. 
We  analyzed  the  performance  of  18,650  cells  as  a  function 
of  number  of  cycles  in  order  to  develop  semi-empirical 
correlations  for  the  state  of  charge  (reversible  lithiation 
capacity)  and  the  battery  resistance  (film  and  polarization 
resistance)  as  a  function  of  number  of  cycles. 

Previously,  based  on  the  performance  data  (charge  and 
discharge  characteristics)  and  the  destructive  physical  anal¬ 
ysis  (DPA)  of  fresh  and  cycled  electrode  materials,  we  iden¬ 
tified  the  main  factors  responsible  for  capacity  fade  of  Li-ion 
systems  with  cycling  [8,9].  According  to  our  analysis,  the 
capacity  fade  can  be  split  into  three  components.  The  first 
part  deals  with  the  loss  in  capacity  due  to  increase  in  resis¬ 
tance  at  both  electrodes.  The  second  part  deals  with  loss  of 
lithiation  capacity  at  both  electrodes.  The  third  part  arises 
due  to  loss  of  active  material  Li+  in  the  cell.  Any  model  that 
will  be  used  to  analyze  the  capacity  fade  has  to  account  for 
capacity  losses  due  to  these  three  terms.  In  this  paper,  we  fol¬ 
low  changes  in  state  of  charge,  film  resistance  and  diffusion 
coefficient  to  quantify  the  capacity  fade.  The  state  of  charge 
accounts  for  the  active  material  loss  (both  primary  and  sec¬ 
ondary  active  material)  while  the  film  resistance  controls  the 
drop  in  the  cell  voltage  with  continuous  cycling.  Since  the 
rate  capability  loss  is  due  to  transport  limitations  induced 
with  cycling  by  formation  of  oxide  layers  on  both  electrodes, 
the  solid  phase  diffusion  coefficient  of  the  limiting  electrode 
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was  considered  as  parameter  to  account  for  the  rate  capabil¬ 
ity  losses.  Discharge  curves  at  different  cycles  were  fitted 
to  first  principles  model.  From  this  fitting,  semi-empirical 
correlations  were  developed  which  account  for  the  varia¬ 
tion  of  the  three  parameters  and  hence  the  battery  capacity 
fade  with  cycling.  These  correlations  were  used  along  with 
a  first  principles  based  dual  Li-ion  intercalation  model  to 
simulate  the  performance  of  the  Li-ion  cell  at  different 
discharge  rates. 

2.  Experimental 

All  experimental  data  that  includes  full-cell  and  half-cell 
studies  for  developing  the  semi-empirical  models  were  ob¬ 
tained  from  our  previous  work  on  capacity  fade  studies  of 
Li-ion  cells  cycled  at  elevated  temperatures  [8,9].  A  brief 
overview  of  the  experimental  work  is  included  here  for 
clarity.  Sony  18,650  cells  with  1.8  Ah  rated  capacity  were 
used  for  all  cycling  studies.  For  charging  the  cells,  the 
conventional  constant  current-constant  voltage  (CC-CV) 
protocol  was  adopted.  A  direct  current  of  1  A  is  used  to 
charge  the  cell  during  the  constant  current  part  and  the  cut 
off  voltage  was  set  to  be  4.2  V.  Subsequently  the  voltage 
was  held  constant  at  4.2  V  till  the  current  drops  to  50  mA. 
Charge-discharge  studies  were  carried  out  in  the  potential 
range  of  2. 0-4. 2  V.  Arbin  BT-2000  battery  cycler  has  been 
used  for  all  cycling  studies.  Three  cycle  numbers  150,  300 
and  800  were  chosen  to  analyze  the  cycled  cells.  Rate  ca¬ 
pability  studies  were  performed  for  the  fresh  cell  as  well  as 
for  the  cycled  cells.  The  cells  were  charged  using  the  same 
CC-CV  protocol  as  used  for  the  cycling  studies  followed 


by  discharging  at  different  rates  (C/9  to  1C).  The  following 
studies  were  performed  on  the  individual  electrodes  of  the 
cell.  The  can  of  cycled  Sony  US18650S  cells  was  carefully 
opened  at  fully  discharged  state  in  a  glove  box  filled  with 
ultra  pure  argon.  Next,  pellet  or  disc  electrodes  were  made 
from  the  positive  and  negative  electrodes  of  Sony  18650 
cell  and  were  used  as  working  electrodes  in  the  T-cell.  Pure 
lithium  metal  was  used  as  the  counter  and  as  reference 
electrode.  Pellet  electrodes  from  the  cells  after  different  cy¬ 
cles  were  cycled  at  very  low  rates  (~C/15)  and  the  intrinsic 
capacity  of  both  positive  and  negative  electrode  was  deter¬ 
mined.  The  difference  in  capacity  between  the  first  cycle  and 
the  present  cycle  yields  the  loss  in  lithiation  capacity  at  both 
electrodes. 

2.7.  The  Li-ion  intercalation  model 

Fig.  1  shows  a  schematic  representation  of  a  typical 
Li-ion  cell  consisting  of  three  regions  namely  negative 
electrode  (graphite),  separator  (poly-propylene)  and  posi¬ 
tive  electrode  (LiCo02).  Both  the  graphite  and  LiCo02  are 
porous  composite  insertion  electrodes.  During  discharge, 
the  lithium  ions  de-intercalate  from  the  negative  electrode 
and  intercalates  into  the  positive  electrode  through  a  highly 
conductive  electrolyte  solution.  As  a  result  of  discharge,  the 
cell  voltage  decreases  since  the  equilibrium  potentials  of  the 
two  electrodes  depend  on  the  lithium  concentration  on  the 
surface  of  the  electrode  particles.  During  charging  of  Li-ion 
cells  the  reverse  reactions  occur.  To  explain  these  phenom¬ 
ena  mathematically,  one  needs  to  consider  model  equations 
that  describe,  mass  transport  of  lithium  in  solid  phases, 
mass  transport  of  lithium  ions  in  the  solution  phase,  charge 


Porous  Composite 
Negative  Electrode 
(Graphite) 

Separator 

Porous  Composite 
Positive  Electrode 
(LiCo02) 

L"  Ln  \ 

^  Ls  ^ 

|  u  q  discharge  >  JJ+  +  g~  +  Q 

Li+  +  e~  +  0  >  Li  0  | 

v  =  0  x=L 


► 

x 


Fig.  1.  Schematic  of  a  typical  Li-ion  cell  sandwich. 
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transport  in  the  solid  phases  and  charge  transport  in  the 
solution  phase  [10,11]. 

The  governing  equations  for  potential  distribution  in  solid 
and  solution  phases  were 

V  •  (CTeffV0j)  -  J  =  0  (5) 

V  •  V02)  +  V  •  (kd  V  In  C2)  +  7  =  0  (6) 


respectively,  where  the  appropriate  definitions  for  effective 
conductivities  (aeff,  /ceff)  and  the  diffusional  conductivity 
(/cd)  were  given  in  detail  in  [10,1 1].  The  model  equation  that 
describes  the  solid  phase  lithium  concentration  is  given  by 


and  for  explaining  the  mass  transport  of  lithium  ions  in  the 
solution  phase  the  following  equation  used  was 

dc2  1  —  t+ 

£23t  =  v'(°2 Vc2)  +  — (8) 

The  equations  describing  the  electrochemical  reactions, 
mass  transport,  and  other  physical  processes  within  the  cell 
are  discussed  in  detail  in  [10-12]. 

Butler- Volmer  kinetic  expression  was  used  to  describe 
the  charge  transfer  processes  occurring  across  both  the 
electrode-electrolyte  interfaces.  Thus,  the  local  volumetric 
transfer  current  density  due  to  charge  transfer  is  given  by 


J  —  cijioj 


exp 


«a 

RT 


exp 


(9) 


where  ioj  is  the  concentration  dependent  equilibrium  ex¬ 
change  current  density  at  an  interface  and  is  given  by 


max 

1 J 


j  =  n,  p  (10) 


The  over-potential  for  the  intercalation  reaction  can  be  ex¬ 
pressed  as 

J 

rij  =  4>1  -  4>2  -  t//,ref - Rf,  j  =  n,p  (11) 

an 


The  equilibrium  potentials  (£//>ref)  of  positive  and  neg¬ 
ative  electrode  are  expressed  in  terms  of  functions  of 
state-of-charge: 


uf  f(0 p) 
Uf  ->  f(6n) 


(12) 


Fig.  2  presents  the  open  circuit  potential  (OCP)  curves  of 
positive  and  negative  electrodes  of  Li-ion  cells  obtained  ex¬ 
perimentally  using  a  T-cell  assembly  with  pure  Li  metal 
foil  as  counter  and  reference  electrode  where  the  electrode 
materials  were  subjected  to  a  very  low  rate  lithiation.  The 
empirical  expressions  (Eqs.  (A.l)  and  (A. 2))  that  represent 
the  variation  of  tAef  as  a  function  of  0 ,  were  obtained  from 

J  ** 

fitting  the  experimental  data.  The  same  empirical  relations 
were  used  for  all  cycling  simulations  of  semi-empirical  mod¬ 
els  presented  in  this  paper. 

The  design  adjustable  parameters  and  other  parameters 
for  the  electrodes  that  are  necessary  for  the  model  are 
given  in  Table  1 .  This  nonlinear  system  of  five  independent 
governing  equations  and  five  dependent  variables  (c\ ,  C2, 
01,  02?  J)  is  solved  as  a  1D-2D  coupled  model  for  the 
three  domains  (negative/separator/positive)  using  FEMLAB 
software. 

In  general,  the  theoretical  specific  capacity  of  any  elec¬ 
trode  material  could  be  calculated  based  on  Faraday’s  law 
that  can  be  expressed  as 

2th,i(As/g)  =  /  =  n,  p  (13) 


Specific  Capacity  (mAh/g) 

0  100  200  300 

4.20 


3.86 

8 

S)  3.52 


> 

=5  3.18 
V 


2.84 


2.50 

0  40  80  120  160 

Specific  Capacity  (mAh/g) 

Fig.  2.  OCP  curves  of  positive  (LiCo02)  and  negative  (graphite)  electrodes  of  Sony  18650  cells. 
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Table  1 


Parameters  for  the  positive  and  negative  electrode  materials 


Symbol 

Units 

Anode 

(graphite) 

Cathode 

(LiCo02) 

U 

[xm 

88 

80 

CTeff 

S/m 

5.764 

12.117 

[xm 

2 

2 

£2 

0.485 

0.385 

max 

C1 

mol/m3 

30555 

51555 

e° 

0.74a 

0.5 

0.83b 

Di 

m2/s 

3.9  x  1CT14 

1.0  X  10"14 

k 

(A/m2 )/ (mol/m3  )3  /2 

1.89  x  10“6a 

9.81  x  10“7 

2.2  x  10“6b 

a  a 

0.5 

0.5 

ac 

0.5 

0.5 

C2 

mol/m3 

1000 

d2 

m2/s 

7.5  x 

10“10 

t+ 

0.363 

«i 

Q  m2 

0.01 

0 

a  Based  on  initial  approach  with  all  capacity  loss  attributed  to  carbon 
electrode. 

b  Based  on  fitting  experimental  half-cell  data  at  low  rates. 


and  the  total  theoretical  capacity  per  unit  area  of  the  elec¬ 
trode  material  can  be  defined  as 

2th,!  (As/m2)  =  F  x  Cmax,!  x  Li,  i  =  n,  p  (14) 

where  Cmax  is  the  maximum  concentration  of  the  electrode 
in  a  phase  which  is  the  property  of  the  material  and  is  given 
by 

Cmax,!  (mol/m3)  =  -P— ,  i  =  n,  p  (15) 

MW/ 

Based  on  the  volume  fraction  of  the  phase,  the  maximum 
capacity  available  can  be  calculated  as 

(2 max,/  (As/m  )  =  F  x  Si  X  Cmax,/  x  Li  (16) 


For  a  completely  charged  Li-ion  cell,  the  negative  electrode 
(graphite  or  coke)  would  be  lithiated  while  the  positive  elec¬ 
trode  (LfiVh^CL  or  LiCo02)  would  be  de-lithiated.  During 
discharge  the  Li-ions  de-intercalate  from  the  carbon  and 
intercalate  into  the  positive  electrode  through  a  highly  con¬ 
ductive  electrolyte  (LiPF6  in  EC/DMC).  Thus,  during  dis¬ 
charge,  the  state  of  charge  of  the  negative  electrode  decreases 
while  the  state  of  charge  of  the  positive  electrode  increases. 
Opposite  will  occur  during  charge.  In  general,  the  SOC  of 
the  negative  and  the  positive  electrode  materials  during  dis¬ 
charge  were  estimated  by  using  the  following  expressions: 


On  =  0°n  - 
6P  =  0p°  + 


^app^ 

Gmax,n 

^app^ 

Gmax,p 


(17) 


where  the  initial  values  of  SOC  of  negative  (0°)  and  positive 
(0p)  electrodes  were  given  by 


Qn 


Thus,  SOC  of  any  electrode  is  defined  as  the  measure 
of  primary  active  material  content  in  it.  For  example,  the 
SOC  of  completely  lithiated  graphite  material  is  100%  and 
it  is  0%  for  completely  de-lithiated  one.  For  developing 
any  capacity  fade  model,  the  parameters  to  be  considered 
for  accounting  capacity  loss  with  cycling  are,  SOC  of  the 
electrode  material,  maximum  concentration  in  a  phase, 
thickness  of  the  electrode  material  and  the  volume  fraction 
of  the  phase  (solid/solution). 

Capacity  fade  correlations  in  literature  attribute  the  entire 
capacity  fade  during  to  the  loss  of  the  active  material.  Since 
Cmax  represents  the  property  of  the  electrode  material  while 
the  thickness  and  volume  fraction  depend  on  the  preparation 
procedures,  it  would  be  a  valid  assumption  to  consider  that 
the  only  variable  that  accounts  for  the  capacity  fade  with  cy¬ 
cling  is  the  SOC.  Since  the  negative  electrode  (carbon)  limits 
the  cell  capacity,  variations  of  lithiation  capacity  of  negative 
electrode  are  considered  for  simulating  capacity  fade. 

Fig.  3  shows  the  discharge  performance  of  Sony  18650 
cells  initially  and  after  500  cycles.  Two  prominent  changes 
observed  in  Fig.  3  are  the  decrease  in  the  overall  cell  capac¬ 
ity  and  the  increase  in  the  voltage  drop  during  discharge.  In 
our  model,  the  SOC  of  the  negative  electrode  is  considered 
for  the  first  change  namely  the  capacity  loss  with  cycling. 
The  second  change,  that  is  the  drop  in  the  cell  voltage,  is  due 
to  the  increase  in  the  cell  resistance.  It  has  been  generally 
accepted  that  a  number  of  unwanted  side  reactions  occur 
during  cycling.  These  reactions  result  in  formation  of  prod¬ 
ucts  that  passivate  the  surface  of  the  electrode.  For  exam¬ 
ple,  pure  Li  metal  has  a  conductivity  of  (1.078  x  107  S/m), 
while  one  of  the  products  (Li2C03,  formed  as  a  result  of 
solvent  reduction  during  charging),  has  a  very  poor  conduc¬ 
tivity  (1.2  x  10  6  S/m).  Thus,  it  is  necessary  to  introduce 
an  additional  term  to  the  over-potential  expression  that  will 
represent  the  potential  drop  at  the  interface  due  to  presence 
of  resistive  film  formed  on  the  electrode  particles.  Thus,  the 
over-potential  can  be  represented  as 

J 

T)  —  01  02  Lref  Rf 

a 


Gth,n 


(18) 


(19) 
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Table  2 


Values  of  6n  and  Rf  for  different  cycle  numbers 


Cycle  number 

< 

Rf  (£2m2) 

1 

0.72 

0.01 

50 

0.672 

0.022 

100 

0.632 

0.0245 

150 

0.62 

0.0271 

300 

0.559 

0.0365 

500 

0.514 

0.044 

The  over-potential  term  given  in  Eq.  (19)  was  used  in 
the  Butler- Volmer  kinetic  expression  (Eq.  (9)).  In  the 
semi-empirical  model  developed  in  this  study,  both  the  SOC 
and  the  film  resistance  of  the  negative  electrode  were  used 
as  an  adjustable  parameter  to  fit  the  experimental  data  for 
various  cycle  numbers. 

The  next  step  involves  development  of  semi-empirical 
correlations  for  two  capacity  fade  parameters  used  to  pre¬ 
dict  the  performance  of  the  battery  as  a  function  of  number 
of  cycles.  The  Li-ion  intercalation  model  was  used  for  all 
simulations.  By  using  two  adjustable  parameters  (the  state 
of  charge  and  the  resistance),  the  discharge  curves  (voltage 
versus  time)  were  simulated  using  the  intercalation  model 
and  were  fitted  with  the  experimental  data.  From  this  fitting, 
the  state  of  charge  (0n)  and  the  film  resistance  (Rf )  were  ex¬ 
tracted  for  different  cycle  numbers.  Table  2  summarizes  the 
values  of  0n  and  Rf  used  to  fit  the  discharge  curves  for  dif¬ 
ferent  cycles.  Based  on  the  variations  of  these  parameter  val¬ 
ues,  the  following  correlations  were  developed  for  the  state 
of  charge  and  resistance  as  a  function  of  the  cycle  number. 

The  decrease  in  the  SOC  of  the  negative  electrode  with 
cycle  number  is  fitted  with  a  one-parameter  correlation  given 
by 

=  <9°  -  k\  (cycle) 1/2  (20) 

where  the  constant  k\  depends  on  cycling  conditions  and 
the  type  of  negative  electrode  material  used.  The  square-root 
dependence  for  the  cycle  number  can  be  accounted  for  by  a 
parabolic  growth  mechanism  for  the  continuous  formation 
of  a  thin-film  solid  electrolyte  interface  (SEI)  layer  over  the 
surface  of  the  negative  electrode. 

The  same  type  of  semi-empirical  correlation  was  found 
to  be  sufficient  to  explain  the  increase  in  the  film  resistance 
with  cycling  and  is  given  by 

R^  =  R°  +  Jt2  (cycle)1/2 


where  k 2  is  a  constant  that  and  R°  is  the  resistance  of  the  SEI 
layer.  The  values  of  the  constants  used  in  the  semi-empirical 
correlations  for  the  film  resistance  and  for  the  SOC  for  two 
temperatures  are  given  in  Table  3.  The  variation  of  the  ad¬ 
justable  parameters  with  cycling  and  the  corresponding  em¬ 
pirical  fit  are  shown  in  Fig.  4.  Based  on  these  correlations, 
the  capacity  loss  and  the  discharge  performance  could  be 
predicted  for  any  cycle  number  by  assuming  that  the  same 
trend  follows  all  through  cycling.  Fig.  5  shows  the  simulated 
discharge  curves  for  several  cycles  with  the  incorporation  of 
the  correlations  developed  for  0n,  Rf.  For  comparison,  the 
experimental  cycle  data  for  cycle  numbers  1  and  800  are 
also  presented. 

Due  to  a  continuous  increase  in  the  film  resistance  dur¬ 
ing  cycling,  the  rate  capability  of  the  cells  decreases.  After 
500  cycles  one  observes  a  clear  difference  between  the  ca¬ 
pacities  obtained  at  very  low  rates  of  discharge  (<C/5)  and 
the  capacity  obtained  with  regular  cycling  rates  (C/2,  1C). 
The  actual  active  material  loss  must  be  calculated  by  using 
only  very  low  discharge  rates.  The  difference  in  the  capac¬ 
ity  obtained  at  low  and  at  high  discharge  rates  will  result  in 
estimation  of  rate  capability  losses.  Since  in  the  above  treat¬ 
ment,  the  entire  capacity  loss  for  C/2  rate  cycling  was  ac¬ 
counted  as  an  active  material  loss  (SOC),  the  model  results 
over  predicted  the  capacity  loss  resulting  from  the  active 
material.  Thus,  it  is  necessary  to  modify  this  model  to  ac¬ 
count  for  the  rate  capability  losses  separate  from  the  active 
material  loss. 

To  further  refine  our  present  approach,  the  capacity  loss 
with  cycling  is  split  to  account  for  both,  the  active  material 
losses  and  the  rate  capability  losses.  The  solid  phase  dif¬ 
fusion  coefficient  of  the  limiting  electrode  is  the  adjusting 
parameter,  which  is  used  to  determine  the  rate  capability 
losses.  The  same  empirical  correlation  used  for  the  variation 
of  film  resistance  with  cycling  is  also  used  here. 

In  our  previous  work  [9],  we  identified  the  following  pa¬ 
rameters  to  be  critical  in  determining  the  capacity  fade  of 
Sony  18650  cells:  the  rate  capability,  secondary  active  ma¬ 
terial  (LiCo02/Carbon)  and  primary  active  material  (Li+) 
losses.  Thus,  the  total  capacity  fade  as  ‘Of  can  be  quantified 
in  terms  of  the  three  factors: 

Q  =  Qi  +  Gn  +  Gin  (22) 

where  Q\  represents  the  rate  capability  loss;  Q\\  the  capacity 
loss  due  to  secondary  active  material  and  Q\\\  refers  to  ca¬ 
pacity  loss  due  to  primary  active  material.  Thus,  the  model 
considers  the  following  parameters:  (i)  the  state  of  charge  of 


Table  3 


Values  of  the  constants  of  the  semi-empirical  correlations  for  him  resistance,  SOC  and  diffusion  coefficient  for  two  different  cycling  temperatures 


Cycling  temperature  (°C) 

Film  resistance  parameter 

SOC  parameters 

Diffusion  coefficient  parameters 

k 2  (£2m2  per  cycle1/2) 

Qi 

O 

3 

k3  (cycle  2) 

k/f  (cycle  2) 

k5  (m2/s) 

ke,  (cycle) 

25 

1.5  x  icr3 

0.837 

8.5  x  icr8 

2.5  x  icr4 

6.134  x  10“17 

1.25  x  103 

50 

1.7  x  10“3 

0.839 

1.6  x  10“6 

2.9  x  10“4 

3.902  x  10“16 

6.91  x  102 
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Fig.  4.  Variation  of  the  adjustable  parameters  6n  and  Rf  with  cycling  and  the  respective  empirical  fits. 


the  limiting  electrode  which  accounts  for  capacity  loss  due 
to  primary  and  secondary  active  materials  (2n  and  2m);  (ii) 
solid  phase  diffusion  coefficient  of  the  limiting  electrode  to 
account  for  capacity  loss  due  to  rate  capability  ( Q\ )  and  (iii) 
film  resistance  to  account  for  the  increase  in  the  polarization 
and  charge  transfer  resistance  with  cycling. 

As  a  result  of  unwanted  side  reactions,  there  is  a  contin¬ 
uous  loss  of  the  active  material  for  every  charge/discharge 
cycle  and  hence  the  kinetic  expression  for  active  material 
losses  could  be  written  as  SOC  varying  with  cycle  number. 
Since  the  negative  electrode  was  limiting  electrode  all  the 
SOC  calculations  were  based  on  graphitic  carbon.  The  ini¬ 
tial  state  of  charge  of  the  negative  electrode  was  estimated  at 
very  low  discharge  rate  (~C/9)  of  a  fresh  cell.  First,  the  spe¬ 
cific  capacity  (2 1  c/9)  was  calculated,  since  the  total  amount 
of  intercalation  host  material  was  known  and  the  SOC  was 


estimated  by  dividing  the  specific  capacity  with  the  theoret¬ 
ical  capacity  of  graphitic  carbon  (372mAh/g): 


(23) 


From  the  active  material  losses  (both  primary  and  secondary) 
found  experimentally  for  several  cycle  numbers,  the  loss  in 
the  SOC  was  calculated  by 

£)lost  2  lost  /O/IN 

9"  =  aT  a4> 


Here  the  term  2lost  includes  both  primary  (2m)  and  sec¬ 
ondary  (2n)  active  material  losses.  Thus  the  SOC  of  the 
limiting  electrode  at  any  cycle  number  is  given  by 

=  e°n  -  elost  (25) 


Fig.  5.  Simulated  discharge  curves  with  the  incorporation  of  empirical  correlations  for  0n  and  R{  to  the  Li-ion  intercalation  model.  The  dotted  curves 
represent  the  experimental  data  obtained  for  Sony  18650  cells  for  cycles  1  and  800. 
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Fig.  6.  Variation  of  SOC  of  negative  electrode  for  the  semi-empirical 
model  with  cycling  at  25  and  50  °C. 


The  capacity  loss  was  quantified  experimentally  for  four  cy¬ 
cle  numbers  (1,  150,  300  and  800)  and  hence  the  correlation 
for  variation  of  SOC  with  cycling  was  developed  based  on 
these  data.  The  rate  of  change  of  state  of  charge  of  electrode 


material  is  expressed  as 

=k3N  +  k4 
dN 

(26) 

initially  0n  =  0°(T) 

(27) 

where  N  is  the  cycle  number.  Fig.  6  presents  the  semi- 
empirical  fit  based  on  Eqs.  (26)  and  (27)  for  the  variation 
of  SOC  of  the  negative  electrode  with  cycling  under  dif¬ 
ferent  temperatures.  The  constant  k 3  accounts  for  capacity 
losses  that  increase  rapidly  under  adverse  conditions  such 
as  cycling  at  high  temperature  and  k\  is  a  factor  to  account 
for  capacity  loss  under  usual  conditions  of  charge  and  dis¬ 
charge.  Table  3  summarizes  the  values  for  the  constants  £3 
and  k\  for  cycling  temperatures  namely  25  and  50  °C.  While 
the  constant  k\  remains  almost  the  same  for  both  cases,  the 
constant  £3,  which  was  stated  as  a  factor  that’s  takes  into 
account  the  cycling  temperature,  increases  with  increase  in 
temperature.  Based  on  the  above  rate  expression  it  is  possi¬ 
ble  to  predict  the  capacity  of  the  cell  at  any  cycle  number  at 
lower  rates  of  discharge.  Still  the  rate  capability  losses  and 
the  polarization  resistance  need  to  be  accounted  to  predict 
the  discharge  curve  at  any  cycle  as  well  as  the  cell  capacity 
at  moderate  to  high  discharge  rates. 

The  performance  studies  on  Sony  18650  cells  [9]  clearly 
showed  that  the  rate  capability  of  the  cells  continues  to 
decrease  with  cycling.  The  fresh  cell  showed  only  a  small 
decrease  in  capacity  with  increase  in  discharge  current  in¬ 
dicating  that  the  cell  possesses  a  very  good  rate  capability. 
Rate  capability  measurements  after  150  and  300  cycles  for 
the  cells  cycled  at  RT  conditions  showed  similar  profiles 
as  that  of  the  fresh  cell.  After  300  cycles  there  is  a  no¬ 
table  difference  in  discharge  capacity  obtained  at  low  and 


high  rates.  For  cycling  at  high  temperatures  (50  °C),  rate 
capability  limitations  set  in  earlier. 

Decrease  in  the  rate  capability  with  cycling  could  be  ex¬ 
plained  based  on  the  film  formation  over  the  electrode  sur¬ 
face  as  a  result  of  the  side  reactions.  Li-ion  diffuses  in  and 
out  of  the  electro  active  material  at  slower  rate  as  a  result  of 
a  continuous  film  formation.  Thus  the  rate  capability  could 
be  treated  as  a  diffusion  limited  problem  and  correlations 
can  be  developed  for  the  change  in  diffusion  coefficient  with 
continuous  cycling.  This  approach  would  be  appropriate  to 
get  better  predictions  of  discharge  performance  at  moderate 
to  high  rates. 

The  diffusion  coefficient  of  negative  electrode  and  the  film 
resistance  were  taken  in  this  model  as  adjustable  parame¬ 
ters.  The  experimentally  obtained  discharge  curves  for  dif¬ 
ferent  cycles  were  fitted  using  the  Li-ion  intercalation  model. 
From  this  fitting  the  value  of  diffusion  coefficient  and  the 
resistance  were  extracted  as  a  function  of  the  cycle  number. 
Next,  using  the  correlation  for  variation  of  SOC  with  cycle 
number,  the  diffusion  coefficient  of  negative  electrode  and 
the  film  resistance  as  adjustable  parameters,  the  discharge 
curves  were  predicted  for  different  cycle  numbers. 

As  stated  earlier,  the  capacity  loss  in  the  model  was  ac¬ 
counted  by  means  of  active  material  losses  (Qu  and  Qm) 
in  terms  of  SOC  and  rate  capability  losses  (gn)  in  terms  of 
diffusion  coefficient.  The  drop  in  the  cell  voltage  with  cy¬ 
cling  was  accounted  by  adjusting  the  film  resistance  term 
that  appears  in  the  Butler- Volmer  kinetic  equation. 

For  cells  cycled  at  RT  conditions,  the  rate  capability  up  to 
300  cycles  almost  remains  the  same  as  that  of  the  fresh  cells. 
There  was  no  significant  difference  in  the  discharge  capac¬ 
ity  obtained  between  very  low  rates  of  discharge  (C/9)  and 
the  discharge  rate  used  for  cycling  (C/2).  The  rate  capabil¬ 
ity  losses  (capacity  at  C/9  —  capacity  at  C/2)  were  estimated 
as  28  and  38  mAh  after  150  and  300  cycles,  respectively, 
which  are  comparable  to  the  rate  capability  loss  of  25  mAh 
for  a  fresh  cell.  Thus  with  the  same  value  of  diffusion  co¬ 
efficient,  we  were  able  to  fit  the  discharge  curves  between 
cycles  1  and  300.  Since  the  rate  capability  continues  to  di¬ 
minish  with  continued  cycling,  the  diffusion  coefficient  had 
to  be  adjusted  or  decreased  to  fit  the  other  discharge  curves 
(cycle  numbers  >  300).  Based  on  these  observations,  the 
correlation  for  variation  of  diffusion  coefficient  with  cycle 
number  was  expressed  as 

dIn  =  ks  exp  (28) 

where  k$  and  k 6  are  constants  that  are  functions  of  both  SOC 
and  cycling  temperature.  Similar  empirical  correlation  was 
developed  for  the  cells  cycled  at  50  °C.  Table  3  summarizes 
the  values  of  the  constants  k$  and  k§  for  the  two  cycling  tem¬ 
peratures.  As  shown  in  Table  3,  the  values  of  both  constants 
increase  with  increase  in  cycling  temperature. 

The  same  correlation  as  used  in  the  first  approach  for 
film  resistance  (Rf)  was  used  in  this  model  also  in  order 
to  simulate  the  drop  in  cell  voltage  with  cycling.  Table  3 
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Fig.  7.  Simulated  discharge  curves  based  on  semi-empirical  model  for  temperatures  (a)  25  °C  and  (b)  50  °C.  The  dotted  curves  represent  the  experimental 
data  obtained  for  Sony  18650  cells  for  the  respective  cycle  numbers. 


summarizes  the  values  of  the  constant  k 2  used  for  the  two 
temperatures. 

Fig.  7a  and  b  presents  the  experimental  discharge  per¬ 
formance  of  Sony  18650  cells  at  various  cycles  discharged 
at  25  and  50  °C,  respectively.  Fig.  7a  and  b  also  show  the 
model  fits  for  comparison.  The  semi-empirical  correlations 
obtained  for  the  change  of  the  state  of  charge,  diffusion  co¬ 
efficient  and  the  resistance  as  a  function  of  cycle  number 
were  used  to  extract  the  parameters  necessary  for  model  pre¬ 
dictions. 

As  shown  in  Fig.  7b,  the  accuracy  of  the  fit  diminishes 
with  an  increase  of  temperature.  This  phenomenon  can  be 
explained  by  taking  into  account  that  for  simulating  the  bat¬ 


tery  performance  at  different  cycles  the  same  empirical  ex¬ 
pression  for  open  circuit  potential  has  been  used.  The  er¬ 
ror  in  fitting  experimental  curves  at  50  °C  indicates  that  the 
property  of  the  intercalation  host  material  changes  with  tem¬ 
perature  that  affects  the  shape  of  the  OCP  curve. 

Fig.  8a  and  b  show  the  capacity  of  Sony  18650  cells  at 
different  discharge  rates.  Discharge  curves  obtained  from 
the  semi-empirical  model  are  compared  to  experimental  data 
for  different  cycle  numbers.  This  analysis  has  been  done 
for  two  different  temperatures  namely,  25  and  50  °C.  For 
all  discharge  rates  between  C/9  and  1C,  a  good  fit  was 
observed  between  the  model  results  and  the  experimental 
data. 
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Fig.  8.  Comparison  of  rate  capability  between  experiment  (dots)  and  model 
(solid  lines)  at  various  cycles  for  cells  cycled  at  (a)  25  °C  and  (b)  50  °C. 


While  the  model  presented  here  can  account  for  various 
capacity  losses,  it  is  still  limited  to  the  system  for  which 
empirical  correlations  were  developed  (Sony  18650  cells  in 
this  case).  Similar  correlations  with  different  parameters  will 
have  to  be  derived  for  other  battery  geometries  and  electrode 
chemistries.  In  order  to  develop  a  generalized  capacity  fade 
model  for  Li-ion  battery,  it  is  necessary  to  use  a  first  princi¬ 
ples  approach  that  takes  into  account  the  basic  mechanism 
for  capacity  fade  common  to  all  Li-ion  systems.  The  model 
should  be  capable  of  explaining  the  capacity  loss  under  cy¬ 
cling  conditions  such  as  charging  to  several  EOCV  and  dis¬ 
charging  to  different  DODs  and  charging  and  discharging  at 
different  rates. 


3.  Conclusions 


fade  was  correlated  to  the  SOC.  The  capacity  loss  due  to 
rate  capability  was  not  considered  separately. 

To  account  for  rate  capability,  a  charge  balance  was  ob¬ 
tained  through  analysis  of  half-cell  data  of  fresh  and  cycled 
cells.  This  model  considered  both  the  active  material  loss 
and  the  rate  capability  losses.  Empirical  correlations  were 
developed  for  the  change  in  SOC  of  negative  electrode  ma¬ 
terial  based  on  experimental  observations.  Solid  phase  dif¬ 
fusion  coefficient  of  the  limiting  electrode  is  the  parameter 
considered  to  account  for  the  rate  capability  losses.  An  em¬ 
pirical  correlation  for  the  variation  of  film  resistance  with 
cycling  was  used  to  explain  the  continuous  increase  in  the 
cell  voltage  drop  with  cycling.  With  this  approach  the  total 
capacity  loss  from  both  active  material  and  rate  capability 
losses  were  determined. 
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Appendix  A 


For  the  graphite  electrode: 


uf  =  0.7222  +  0.13876>n  +  0.0296>n/2  - 

0.0019  m  or,  , 

+  — —  +0.2808e(0-90_15e") 
a  1.5 
un 

—0.7984 e(O-44650n-°-41O8) 


0.0172 


(A.l) 


For  the  FiCo02  electrode: 


-4.656  +  88.6696*2  _  4oi.n96»4  +  342.9096*6 
—462.4716**  +433.434 ef 

-1  +  18.9336*2  -79.53204  +  37.31 16*| 
—73.0836**  +  95.960p° 


(A.2) 


The  solution  phase  conductivity  as  a  function  of  concentra¬ 
tion  C2  (in  mol/dm3)  is  [13]: 

Keff  =  k4°  =  (4.1253  X  10"4  +  5.007C2  -  4.7212 

x  lO3^  +  1.5094  X  lO6^  -  1.6018  x  10*4)4° 

(A. 3) 


A  semi-empirical  capacity  fade  model  has  been  developed 
for  Fi-ion  cells.  As  a  first  step,  empirical  correlations  were 
developed  for  state  of  charge  and  film  resistance  based  on 
fitting  the  experimental  data  with  Fi-ion  intercalation  model. 
The  entire  capacity  loss  with  cycling  was  assumed  to  be  as  a 
result  of  active  material  degradation  and  hence  the  capacity 
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